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ABSTRACT 

^ ' We re-examine the constraints which can be robustly obtained from the 

Q I observed temperature function of X-ray cluster of galaxies. The cluster mass 

function has been thoroughly studied in simulations and analytically, but a 
direct simulation of the temperature function is presented here for the first time. 
Adaptive hydrodynamic simulations using the cosmological Moving Mesh Hydro 
^ I code of Pen (1997a) are used to calibrate the temperature function for different 

I popular cosmologies. Applying the new normalizations to the present-day 

' cluster abundances, we find erg = 0.53 iO.OSOo*^'^^ for a hyperbolic universe, and 

CTg = 0.53 ± 0.05fig for a spatially flat universe with a cosmological constant. 
^ I The simulations followed the gravitational shock heating of the gas and dark 

^ I matter, and used a crude model for potential energy injection by supernova 

heating. The error bars are dominated by uncertainties in the heating/cooling 
^ models. We present fitting formulae for the mass-temperature conversions and 

in ' cluster abundances based on these simulations. 



X 



1. Introduction 

X-ray studies of clusters of galaxies have provided a host of quantitative data for the 
study of cosmology. They are ideally suited for the normalization of structure formation 
since their mass (~ IO^^Mq) is close to the mass formed from the collapse of an 8h~^ 
Mpc radius sphere of matter. They provide a direct dynamical measure as, the linearly 
extrapolated variation in mass in such spheres. We can numerically compute their formation 
directly through simulations from linear perturbation theory into the non-linear regime 
with appropriate boundary conditions. Only about a percent of the mass of the universe 
is in virialized rich clusters. For a Gaussian random field the variance depends only 
logarithmically on the abundance, making an accurate measurement possible even in the 
presence of substantial Poissonian errors. 
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While the gravitational mass in bound objects is easily computed using analytic 
approaches such as the Press- Schechter formalism (Press and Schechter 1974) (hereafter 
P-S) or N-body simulations, it is not so easily measured. The most robust measurements 
of the cluster abundance currently rely on the cluster temperature function, the number 
density of clusters above some temperature A^(> kT) expressed in units of /i"^ Mpc^"^. Henry 
and Arnaud (1991) (hereafter referred to as HA) and Edge et al. (1990) have analysed 
the HEAO-1 A2 all sky X-ray survey to obtain temperatures for clusters at a flux limit 
of F2_iokeV > 3 X lO"^-*^ erg/cm^/sec. The sample contains the 25 X-ray brightest clusters 
and is over 90% complete at galactic latitude |6| > 20". By virtue of their brightness 
the temperatures are easily and well determined. In this paper we will make use of this 
excellent objective sample. 

Prom the theoretical stand point, the cluster temperature function is a clean 
measurement since the normalization of fluctuations can be established independently of 
the baryon fraction Vlb, or the Hubble constant Hq = 100/ikm/sec/Mpc. Clusters are close 
to isothermal, both observationally and in simulations, which makes their temperature 
determination robust and insensitive to numerical resolution or telescope angular resolution. 
The temperature of a cluster depends primarily on the depth of the potential well of the 
dark matter, and the state of equilibrium of the gas. This is in contrast to the cluster 
luminosity, which depends strongly on small scale parameters like clumping and core radius. 
While heating or cooling processes could change the temperatures slightly, the gravitational 
potential provided by the dark matter buffers the system and causes the temperatures in 
hydrostatic equilibrium to only be weakly affected by energy injection or loss. We will 
verify this effect through direct simulations of heat injection. Furthermore, measurements 
of the galaxy velocity dispersion appears to be in good agreement with the temperatures, 
with the ratio of gas temperature to velocity dispersion /3fit = kT/fj.nipa'^ = 0.95 ± 0.05 
(Bahcall and Lubin 1994). This suggests that non-thermal pressures are negligible. 

Recent work by Eke, Cole and Frenk (1995) (hereafter ECF), and by Viana and Liddle 
(1996) (hereafter VL) used the cluster temperature function to normalise the amplitude 
of fluctuations in spheres of radius 8h^^ Mpc, denoted a^. They compared the observed 
cluster temperature function of HA to Press- Schechter calculations calibrated by N-body 
simulations. These studies did not include any gasdynamic effects, leaving open questions 
about the statistical properties of cluster equilibria. They concluded that the present 
day normahzation in a flat standard CDM is as — 0.5. Modelling the predicted cluster 
abundance in a subcritical density universe is more difficult. VL propose erg oc Jlo°'^ 
while ECF have a smaller error budget in the exponent. ECF only measured the mass 
function in simulations. In order to convert from a mass function to a temperature function, 
these models all assumed that all clusters are perfect singular isothermal spheres which 
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formed at some prescribed redshift. While the errors in the normahzation dg from the 
Poissonian scatter are very small, any error in mass-temperature conversion translates into 
a similar sized error in the normalization. That issue is the point of largest uncertainty in 
these semi-analytic calculations. In order to address this problem, we extend their work 
by performing full hydrodynamic simulations. We directly measure the X-ray emission 
weighted temperature in the simulation. The simulations are used to directly measure the 
mass-temperature relation which are then fed into the P-S calculations. This approach 
allows us to explore a large dynamic range in volumes and parameters. 

In this paper we will not use the term "open universe", but instead call universes 
with negative spatial curvature "hyperbolic universe". This circumvents the misleading 
suggestion that hyperbolic universes should be spatially infinite (Pen and Spergel 1995). 
We will proceed in section ^ to describe the simulations. In section ^ we summarize the 
observed cluster temperature function. Our new results are presented in section ^. We 
discuss the cosmological implication of these results in section |^. 

2. Simulations 

We use the Moving Mesh Hydrodynamic (MMH) code (Pen 1997a). It implements a 
Total Variation Diminishing (TVD) (Xin and Jin 1995) high resolution shock capturing 
hydrodynamics scheme on an adaptively deforming mesh. The gravitational potential is 
solved using a multigrid iteration (Pen 1995). Dark matter is modeled using a particle- mesh 
algorithm on the same moving grid. The grid is continuously adjusted to maintain an 
approximately constant mass per cell. This is achieved through a pure potential flow 
grid velocity field. The full Euler fluid equations are solved on this moving mesh, and 
the fluid is allowed to develop vorticities. By following the mass, the MMH code has 
improved spatial resolution in dense regions such as clusters of galaxies. Since their cores 
are 10'^ — 10^ times overdense, any mass based method will have a ten to twenty fold better 
length resolution in these high density regions. A recent study has compared this code to 
several other cosmological hydrodynamic codes for a cluster formation scenario (Frenk et 
al. 1997), and good general agreement exists between this code and the others compared. 
Excellent agreement was found for the temperature properties of simulated clusters. Similar 
agreement was obtained in a comparison to a suite of existing codes (Pen 1997a, Kang et 
al. 1994). All simulations were run on an SGI power challenge at the National Center 
for Supercomputing Applications. They all used 128^ grid cells and 256^ particles. The 
initial power spectrum was taken from Bardeen et al. (1986) (hereafter called BBKS). 
The simulations were started at a redshift z = 100. Table |l| summarizes the cosmological 
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parameters used in each simulation. Cooling has not been incorporated in these simulations, 
and will be addressed in a future paper (Cen et al. 1998). 

For the CDM model we used the best fit values suggested by ECF with 
Qq = 1, Qb = 0.05, h = 0.5, as = 0.5, n = 1 on a box of side length 80h~^ Mpc. 
n denotes the unprocessed power spectrum (BBKS), whose Harrison- Zeldovich- Peebles 
(hereafter HZP) value is n = 1. The simulations were repeated two times with different 
random seeds to improve the statistical measures. These comprise models CDMl-3. We 
note that the finite box sizes truncates the long wave modes, which has the effect of lowering 
the fluctuation amplitude in the numerical realization by 15% (Gelb and Bertschinger 1994, 
Pen 1997b). To check its significance, and to test the resolution dependence, we ran a 
simulation in a larger box of side length 200h~^ Mpc which we call CDM4. At such small 
values of Qb the corrections to the power spectrum are negligible compared to other sources 
of error (Holtzman 1989). As long as the gravitational potential is dominated by the dark 
matter, we can rescale the result to any value of Qb easily. We then ran a model with 
heat input from early star formation, which we call PREHEAT. At a redshift 2; = 1, we 
injected 1 keV of energy per nucleon into the plasma. Details are described in section ^. We 
next ran a cosmological constant model Q/<^ = 0.63, erg = 1 in a box of side length 120h~^ 
Mpc. We repeated the simulation with identical initial conditions for a hyperbolic universe 
fio = 0.37. These are called LCDM and OCDM respectively. 

Since clusters are rare objects, their identification is a relatively simple matter. We 
use the gas densities to identify clusters. We search the volume for all density peaks in the 
unsmoothed gas field which are overdense by at least 200 over the mean cosmic density, and 
separated by at least 2h~^ Mpc. If two peaks are closer than that distance, the less dense 
one is discarded. From this sample, we computed the total emission weighted temperature 
for each cluster in a lh~^ Mpc radius, weighting the temperature in each cell by p^T^/^, 
and sorted these clusters by temperature to find the cumulative temperature function. 

3. Observed Temperature Function 

We use the cluster sample from HA. The sample is based on the objects from the 
HEAO-1 A2 survey, which identified all objects at a flux limit > 3 x 10^^^ erg/cm^/sec 
in the 2 — 10 keV band at a galactic latitude \b\ > 20°. The 2-10 keV window over which 
HEAO is sensitive is well matched to the temperatures of rich clusters. HA identified 25 
clusters with that sample. They range in temperature from 2.5 — 9.5 keV, with luminosities 
in the range 10^^ — 10^^ erg/sec. The furthest cluster is at redshift z = 0.09, which 
is sufficiently low that evolutionary corrections are small. This sample has excellent 
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completeness properties and most clusters have multiply measured temperatures. The 
typical temperature error is 10%. Following ECF, we measure the cumulative cluster 
density 

N{>kT)= y: ^ (1) 

kTi>kT "^max,! 

where Knax,i is the maximal volume to which each cluster could have been seen at the flux 
limit. The formal Poisson error on the estimate is given by 

a{>kTr= 77^- (2) 

kTi>kT ''^max.i 

This tends to underestimate the true error for several reasons. Since we are measuring 
a cumulative abundance, errors at different temperatures are correlated. Clusters tend 
to cluster strongly with each other with correlation lengths of about 20/i~^ Mpc (Bahcall 
1996). This boosts the error in (0) for the low luminosity clusters where Kaax is comparable 
to the correlation volume. The errors are Poissonian with respect to the true underlying 
distribution, not with repect to the estimated density. Using the estimated density will 
again systematically underestimate the errors. Instead of attempting to model these 
uncertain errors directly, we will use a heuristic flt based on the estimates of ECF. They 
found that the abundance is best represented near 5 keV, which lies in the center of the 
temperature window. This value is close to the temperature formed by the collapse of an 
8/;,"^ Mpc radius sphere in a = 1 model, and empirically appears to be the pivot point 
of the distribution as one varies parameters. Normalizing abundances at that temperature, 
ECF found statistical errors result in variations of cts of only about 2%. As we will see 
below, the errors arising from thermal history uncertainties and numerical limitations 
are signiflcantly larger, so we will neglect the Poissonian errors hereafter. The cluster 
abundance and various past flts are shown in Figure |I] renormalized using the results of this 
paper (see equations H, ^ and |Tl] below). 



4. Mass- Temperature Relation 

To obtain an analytic estimate of the temperature function, we will use the 
Press-Schechter Ansatz. The fraction of mass in bound objects is 




/(>M) = ^/- l^e--''du. (3) 

<7(M) 



cr(M) is deflned as the RMS density fluctuations in tophat spheres of mass M. The 
distribution function has been multiplied by 2 such that all the mass is accounted for when 
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a oo. 5c = 1.686 is the linearly extrapolated overdensity at which an object virializes. 
Formally it is only exact for Qq = 1, but it varies only by a few percent to Qq = 0.3 and 
we will consider it to be constant. We will numerically normalize the mass-temperature 
relation, which will absorb any changes in 6c and simplifies calculations. We define the 
dimensionless mass m = M/Ms where Mg = 47rp(8/i-^Mpc)V3 is the mass contained in an 
8/i-^Mpc sphere, p = SnoH^/SirG is the mean density of the universe today in terms of 
the Hubble constant Hq = lOO/ikm/sec/Mpc. Differentiating (0) we obtain the differential 
number density of objects dn/dm = {p/Ms'm)df /dm as 

dn 2 p 6c 
dm V Mg am^ 

The observed abundance (|I|) is a cumulative statistic, so it is desirable to integrate the 
differential density function (^). To simplify the algebra, we will assume a pure powerlaw 
dependence a = agm^". a is related to the effective power spectrum P{k) = For 
CDM-like power spectra, the BBKS fit to the power spectrum at r = 8h~^ Mpc has 

a = 0.222 + 0.2495r - 0.0232(^0, (5) 

where F = Qohexpl—Q^ — (^fe/^o)] (VL). (|^) is defined by requiring the ratio of oq/ctx 
to be exact, where X is the radius which forms a 5 keV cluster in the three cosmological 
models which were simulated. We can then integrate (Q) 

Q / X2 \ ^ O- — 1 6'^ 

^ ' 20487r3/2 \^2(T|y ^ 2a 2af ^ ^ ' 

where F[a,a;] is the incomplete F function. The cumulative number abundance (P) has units 
of clusters per Mpc^//i^. The asymptotic semi- convergent sequence for F (Arfken 1985) 
allows us to expand (||) 

n{> m) ^ 2.63 x 10-(//-)^^^ f 1 - ^ + - + + +.■] /.^Mpc^ 

(7) 

where we have abbreviated p = (5^/2(j|, v = pm^^, p = (a + l)/2a. The terms in (0) 
initially converge for large f , but at some point diverge again. Each term has an alternating 
sign, and each approximation brackets the true solution. For our purposes, we can truncate 
the sum and modify the first two terms as 

n(> m) ~ 2.63 x 10-^'/'"^^^^!^ (l - £) /i^^Mpc^l 

The accuracy of is shown in Figure 0, where we have used the mass-temperature relation 
derived below. Three sets of lines are shown, corresponding to an integral of (^) using 
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the BBKS power spectrum (solid lines), the abundance using the power law fit with 
and (|^) (dashed lines) and the two term asymptotic expansion (0) (dotted lines). For this 
illustration, we picked one model with = 1, and one cosmological constant model with 
= 0.65. We see that the approximations are accurate to about 3% in temperature over 
the temperature range 3-10 keV. 

For the next step, we need to convert the cumulative mass function into a temperature 
function. Various approaches, for example ECF, assumed all objects to be perfect isothermal 
spheres forming at ^ = 0. In the chain of calculations, this step is the least certain link, and 
it is the purpose this paper to quantify it. We will continue to rely on the virial relation 
m (X T^/^, which follows if the formation redshift and state of equilibrium is uncorrelated 
with cluster mass. Bryan and Norman (1997) recently demonstrated this relationship to 
hold to a good approximation in hydrodynamic simulations. We then only need to know 
the relation coefficient Tg 

where Tg is the effective temperature of a cluster which forms during the collapse of an 8h~^ 
Mpc radius sphere. In a low density universe the same sphere contains much less mass, and 
consequently Tg depends on Qq. While the overdensity parameter Sc in the top hat model 
has a weak Qq dependence, we will fix its value at Sc = 1.686 and absorb all modeling into 
Tg. 

To solve for Tg we identified clusters in each l^o = 1 simulation, resulting in the 
numerical cumulative distribution n{> T) shown as the step lines in Figure |^. We then 
applied (^ to the linearly extrapolated density field of the simulation's initial conditions 
6p{z = 0) = 6p{zi) X (1 + Zi). This density field is then smoothed on tophat spheres 
of varying radii, from which we obtain the predicted Press-Schechter abundance of the 
collapsed mass fraction shown as the dashed diagonal lines in Figure |[ We then solve for 
Tg such that (H) agrees with the predicted abundance in the simulation volume for the 5 
hottest clusters. This approach compensates for the finite box size effect which results in a 
loss of power on large scales (Pen 1997b). The infinite volume average is the solid diagonal 
curve in Figure |^. We see that the smaller boxes systematically underestimate the cluster 
abundance due to the suppression of cxg from the truncation of the power spectrum in the 
finite box. 

The best resulting fit is 

Tg = 4.9 ± .2keV. (10) 

This value should be compared to Tg = 5.5 keV for the perfect isothermal sphere model 
(ECF). We would expect the ECF model to overestimate Tg for several reasons. The gas is 



- 8 - 



almost certainly not in perfect hydrostatic equilibrium, which will lower the termperature. 
A scatter in the mass-temperature relation also creates more hot clusters due to Malmquist 
bias. Cluster profiles may have departures from isothermality, with slight temperature 
gradients throughout the cluster. The X-ray emission weighted temperature can be slightly 
different from the mean mass weighted virial temperature. ECF provide for a fudge factor 
P which describes the ratio between the temperature they used and the actual statistical 
temperature of clusters. Note that the mean mass-temperature relation for an ensemble of 
clusters is not sufficient to substitute correctly into the Press-Schechter relation (H). Since 
clusters are very rare, the mass function is steep, and any scatter in the mass temperature 
relation will introduce a bias in the abundance. What we really want is to normalize 
directly to the temperature function measured in simulations. One can rescale the ECF 
result by choosing (3 = 1.1 in their temperature conversions. Concurrent work by Bryan and 
Norman (1997) produced results very similar to this study. They defined a parameter Jt as 
the ratio of the actual temperature of a cluster compared to that obtained from the tophat 
model with some assumed radial density profile. In our notation, a rough correspondence 
would be fx = Tg/S.S, for which we would obtain fx = 0.89 ± 0.03, while they found values 
in the range 0.75 ^ fr ^ 0.92. 

The error interval is the 1 — a standard deviation comparing the different simulations 
and represents departures from a deterministic P-S theory. To first order, we have 
compensated for numerical cosmic variance, and also the loss of large scale power from the 
finite box size. The simulations only have a limited sampling volume, and do not directly 
simulate the rarest, richest clusters. By using the P-S formalism, we can extrapolate the 
simulation normalization to larger volumes by assuming the mass temperature relation 
scales as one would expect from virial equilibrium. Each simulated model still had at least 
one cluster above the pivot temperature of 5 keV. In this analysis, P-S allows us to reduce 
the error bars by simultaneously fitting a range in cluster temperatures. 

The first non-gravitational effect which needs to be incorporated is the effect of 
heat injection from stars. Direct observations of the iron line emission suggests that the 
intra-cluster medium (ICM) is not pristine, and has passed through an earlier generation 
of stars. This may have raised the initial entropy of the gas. The present day metallicity 
of the ICM is near 1/3 solar, from which we may infer up to 1 keV of energy per nucleon 
to have been injected (Loewenstein and Mushotzky 1996). While it is not known when this 
might have happened, clusters at redshift 0.3-1 appear to have similar metalicities as nearby 
clusters (Mushotzky and Loewenstein 1997, Hattori et al. 1997). Early enrichment would 
have a smaller effect since the adiabatic expansion of the universe cools the preheated gas. 
In model PREHEAT we inject gas consistent with the observed lack of evolution to 2; = 1. 
The most extreme model postulates 1 keV of energy injection aX z = 1. We simulated 
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such a model by evolving a simulation to z = 1, raising the thermal energy everywhere 
by 1 keV, and continuing the evolution to 2; = 0. We then measured the present cluster 
temperature function. We find that clusters are slightly hotter with Tg = 5.3 keV. We 
should consider this an upper limit on the plausible effect of heating. Because the gas 
remains in hydrostatic equilibrium with the dark matter potential, the injection of 1 keV 
only raised the mean temperatures by 0.4 keV. For our choice of flat universe parameters, 
this lowers the normalization to erg = 0.50 relative to the adiabatic value of 0.53 (see below). 
Future work will also address the effects of cooling (Cen et al. 1998). Fabian (1994) 
estimates that up to 20% of the X-ray luminosity of a cluster may arise from a cooling flow, 
which might affect the emission weighted temperature by a similar amount as the heating. 
The latter depends sensitively on the smallest scale inhomogeneities, and poses a larger 
computational challenge. For now we will assume temperature errors to be symmetric. 

Once we have solved for Tg using the fio = 1 simulations, we proceed to repeat the 
procedure for the A and hyperbolic cosmologies. By using the same random seed, we 
expect the differences between the models to be modeled more accurately than each model 
individually. The VL dependence is incorporated into to mass-temperature relation as 



The 2/3 scaling accounts for the smaller virial mass enclosed in the '^h^^ Mpc spheres as the 
density is lowered. 7 will be used to parameterize remaining corrections, such as the change 
in formation redshift, change in virial radius and collapse density. Solving for 7 from the A 
simulation yields 7 = 0.283 for a A universe and 7 = 0.133 in the hyperbolic scenario. The 
temperature function flts for the simulated low density parameters are shown in Figure 

We now have all the required relations to solve for ag given the cluster abundance. For 
a given VIq, we have the left hand side of Equation (^ from HA at any given temperature, 
for which we follow ECF and use kT = 5 keV. Equations (^ and (|TTp allow us to convert the 
temperature into m. The only remaining unknown variable is (Xg which determines ^ and v. 
We then obtain the result ag = 0.53 ± 0.05^0 in the A model, and ag = 0.53 ± 0.051^0 
in the hyperbolic universe. The error bars are obtained by linearly adding the uncertainty 
in the mass-temperature relation to the effects of supernova heat injection. Figure |^ 
summarizes the results from of the numerical normalizations. In the hyperbolic model, 
matter domination ends at a relatively higher redshift 2; ~ (l/fi) — 1 than in a cosmological 
constant model. Cluster of the same mass will thus have a smaller radius in the hyperbolic 
case, and thus a higher virial temperature. This accounts for the smaller value of 7 and (Xg 
in hyperbolic models. Our modeling of the temperature function ([n|) implicitly accounts 
for all these effects since it is normalized by simulations. 




(11) 
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5. Cosmological Implications 

In the context of structure formation, the simplest adiabatic models, which could arise 
for example from inflation, have only a single free parameter to normalize the spectrum 
of fluctuations for a HZP spectrum. For a COBE normalized flat universe, this implies 
(78 = 1.2, which is at great odds with the observed cluster abundance. If we wish to retain 
a flat universe, one can can try to change one of several parameters. The first would be 
the baryon fraction (White et al. 1995) which suppresses fiuctuations due to acoustic 
oscillations and Silk damping. COBE normalized CDM agrees with the cluster abundance 
for Vti) = 0.45, which requires a dramatic revision of Big-Bang nucleosynthesis (Walker et al. 
1991). A different parameter which can lower the COBE normalized value of ag is a tilted 
power spectrum, i.e. deviations from HZP. Leaving only this one parameter free, we flnd a 
satisfactory flt for n = 0.63 when we also allow for tensor modes. This violates the limits 
on the slope allowed by the 4 year COBE data (Wright et al. 1996). A combination of 
these two parameters violates each of these constraints more weakly, and one could envision 
combinations such ets n — .7, = .15, which would also be consistent with the cluster 
gas fractions (White et al. 1993, White and Fabian 1995). One can also take more radical 
departures, and lower the Hubble constant to ~ 0.3 (Bartlett et al. 1995). We conclude 
that no single parameter modification of O = 1 infiationary cosmology is even marginally 
consistent with observations. 

For low flo models, the opposite problem arises. COBE normalized fiuctuations result 
in too low values of ag. For the hyperbolic universe with parameters Qq = 0.37, h = 0.7, 
we obtain cTg = 0.57, significantly lower than the 0.83 suggested by the cluster abundance. 
To raise it to match the cluster abundance, one can either raise Qq to 0.43, raise the 
Hubble constant to /i = 1.1, or introduce a tilt n = 1.17, or any combination thereof. The 
cosmological constant models also have their share of free parameters. By lowering flo, we 
lower the normalization, but it increases with larger Hubble constant h. For the choices 
under discussion — 0.37, h — 0.7, we have COBE normalized erg = 1.0, just shghtly 
higher than the cluster abundance suggests. This can easily be addressed by lowering h to 
0.63, or using the shght tilt n = .95 as suggested by Ostriker and Steinhardt (1995). 

A third alternative is to consider non-adiabatic initial conditions, for example from 
topological defects (Pen et al. 1997). In these models the P-S abundance must be modified 
to account for the non-Gaussianity (Chiu et al. 1997). These models have not been studied 
directly with hydrodynamical simulations of the cluster temperature function, but the 
preliminary results indicate a cluster abundance consistent with a COBE normalized HZP 
spectrum. We note that these models have problems with other observations, including the 
galaxy power spectrum and small scale microwave background anisotropics. 
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The HA sample has the great advantage that it has well established completeness 
criteria, which allows us to accurately measure cluster abundances. Many more clusters 
have measured temperatures, and one can ask how those might affect our estimates. Most 
importantly, we at not limited by Poissonian statistics but rather by systematic errors. 
It has been proposed that the evolution of the temperature function is a strong test of 
cosmologies. At fixed temperature, the difference in cluster abundance at z = 0.5 is over an 
order of magnitude between flat and hyperbolic models (ECF). Again, we are not limited 
by statistics. Instead, we must ask what the expected difference in cluster temperature 
is at fixed abundance. Presently, the hottest cluster in a 300/i^^ Mpc radius is about 8 
keV (see Figure |l|). At a redshift z = 0.3, the hottest cluster in the same sample volume 
would be essentially the same temperature in a hyperbolic universe, and about 6 keV in a 
flat universe. One must make sure that systematic temperatures measurement errors are 
less than 25% at these redshifts. Since the metallicities have not evolved significantly, it 
appears that heating will not contribute significantly to the cluster temperature evolution. 
The presence of hot clusters such as A2163 and MG2016 (Hattori et al. 1997) may well 
pose problems for flat cosmologies (Pen, David and Tucker 1997). Recently, Carlberg et 
al. (1997) studied the abundance of galaxy clusters at intermediate redshift using galaxy 
velocity dispersions (see also Fan et al. 1997). The difficulty here is converting local cluster 
temperatures into velocity dispersions (Bird et al. 1995, Bahcall and Lubin 1995) since the 
expected change is small. Any scatter in the velocity-temperature relation will introduce 
systematic biases which need to be understood. A systematic overestimate of 10% in the 
line-of-sight velocity dispersion in the CNOC sample relative of the converted (T„ from 
the HA temperature data is sufficient to offset the decrease in the cluster abundance 
predicted in an f2 = 1 cosmological model. It is essential in these comparisons that 
homogeneous samples are used at both low and high redshift which are checked directly 
against simulations. 

Clusters of galaxies provide, coincident ally, a similar constraint on VLq as velocity field 
measurements. Measurements of peculiar velocities constrain j3v = ^o/b, where the bias b 
is the ratio of fluctuation in galaxies relative to the dark matter (Strauss and Willick 1995). 
Typical values for (3^ are in the range 0.3 — 0.8. Cluster abundances from the temperature 
function constrain a very similar function ag^Q^'^ ~ 0.5, where as = 1/b for optical galaxies. 
While velocity fields measure fluctuations in the linear regime, the cluster abundances are 
determined for highly non-linear bound objects. It is reassuring that the values obtained 
from the two very different methods are consistent with each other. The downside is that 
we cannot determine Qq and b independently using present day measurements alone. 



- 12 - 



6. Conclusions 

The most robust cosmological constraints using clusters of galaxies come from the 
cluster temperature function, which depends primarily on the gravitational potential 
wells of the dark matter. We used new gas dynamic simulations to test the N-body and 
Press-Schechter estimates by ECF and VL. We found good general agreement in the 
normalization to the observed cluster temperature function with erg = 0.53 ± 0.05Qq^'^^ 
for a hyperbolic universe, and (Tg = 0.53 ± 0.05fio "'^^ for a spatially fiat universe with a 
cosmological constant. This result is only weakly sensitive to models of the thermal history 
of the intra-cluster medium, which we have modeled with preheating. We have presented 



improved Press-Schechter fits to predict the mass-temperature-relation scalings (11) more 
accurately for a range in values of Qq- Applications to high redshift X-ray clusters are in 
progress (Pen, David and Tucker 1997). 

COBE normalization with a Harrison-Zeldovich spectrum over-predicts cluster 
abundances by many orders of magnitude, but this can be addressed if the spectrum 
is strongly tilted and a large baryon fraction is invoked. The Qq = 0.35 hyperbolic or 
cosmological constant models lie closer to observations on all measures, including the age 
of the universe, the slope of the temperature function, the gas fraction and the slope of the 
galaxy power spectrum just to mention a few. While no single measurement is at a very 
high significance, the combination does appear to carry a heavy vote. Lensing statistics 
(Kochanek 1996) and deceleration parameter measurements (Perlmutter et al. 1997, Pen 
1997c) would favor a hyperbolic universe over one with a cosmological constant, and 
alternative scenarios (for example the string dominated model, Spergel and Pen 1996) are 
also viable. COBE normalized topological defect models fare reasonably well on the cluster 
abundance (Chiu et al. 1997). 
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Cluster Abundance Data 




T (keV) 

Fig. 1. — The stepped solid line is the observed cumulative temperature function from HA. 
The upper and lower thin stepped lines indicate 1 — a Poisson errors. The diagonal solid 
curve represents the cluster abundance normalized to erg = 0.5 (ECF), while the dotted curve 
is the VL normalization (cxg = 0.6). 
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Exact vs Approximate Press-Schechter 




kT (keV) 



Fig. 2. — Comparison between the temperature function obtained by integrating a full power 
spectrum compared to the scale-free power law and its approximation. We see that Equation 
P a good approximation on the scales of interest. 
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Fig. 3. — The cumulative CDM = 0.5 simulation temperature function. The long diagonal 
curve is the Press-Schechter cluster abundance with the best fit Tg = 4.9keV. Tg is the 
temperature of a cluster formed from the collapse of an 8h~^ Mpc radius sphere. The 
stepped lines were obtained from four CDM simulations. The upper left steps are from the 
three simulations CDMl-3 with box size 80h~^ Mpc, and lower right steps are from CDM4 
simulated in a 200h~^ Mpc box. The dashed diagonal lines are the Press-Schechter cluster 
abundance obtained obtained with the actual simulation realizations. For the smaller boxes, 
a loss of power due to the absence of long waves causes the long dashed line to be lower than 
the solid ensemble average. 
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Low Q Models 




T (keV) 

Fig. 4. — The simulation fits to iiypcrbolic and cosmological constant dominated universes. 
Both simulations had identical random seeds and = 0.37 and erg = 1. The dashed 
abundance is for the cosmological constant model. The sohd line is the hyperbolic model. 
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Models Comparison 




3 4 5 6 7 8 9 

T (keV) 



5. — Comparison of the best fit normalizations to the Henry and Arnaud cluster sample. 
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model 


^^0 






h 


L(/i-i Mpc) 


CDMl-3 


1 





0.5 


0.5 


80 


CDM4 


1 





0.5 


0.5 


200 


PREHEAT 


1 





0.5 


0.5 


80 


OCDM 


0.37 





1 


0.7 


120 


LCDM 


0.37 


0.63 


1 


0.7 


120 



Table 1: Simulations used for this study. All simulations used a 128^ grid with 256^ particles. 
All models have a baryonic fraction Qft/i^ = 0.0125 and compression limiter Cmax = 10. The 
effective resolution is (128cinax)^ (Pen 1997a). The two low VL models used identical random 
seeds. 



